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The description of self-gravitating stellar systems as 
equilibrized objects in the microcanonical thermodynam- 
ics applied also to "Small" systems 1] with sharp total 
energy and angular momentum demands a spherical con- 
tainer centered in the center of mass of the total system. 
Of course in reality there is no such container. The con- 
tainer is necessary to avoid particle loss by evaporation. 
The statistical equilibrium thus obtained is a metastable 
equilibrium which lives over times such that the mass loss 
due to evaporation is negligible, and makes sense as long 
as nuclear hydrogen burning stabilizes the stars against 
collapse. 

In order for the calculated equilibrium configuration 
to have anything in common with the real free stellar 
system the effect of the container on the obtained mass 
distribution must be minimized. I.e., only configurations 
with the center of mass fixed at the center of the con- 
tainer are admitted. Moreover, the density near the con- 
tainer walls should be small. This was achieved in our 
calculation Q, H, |^ |^ by setting the contributions of the 
spherical harmonics J bi^i,„i{x)x^dx = 0. 

In the molecular-dynamics simulation by Ispolatov 


and Karttunen this necessary constraint was forgotten. 
Therefore the spatial distribution they show in their Fig- 
ure 1 is irrelevant for reality because the container has a 
non-negligible momentum and the center of mass of the 
gas is not at rest. Certainly after removing this crucial 
restriction it is impossible to find any realistic configu- 
ration like binary stars rotating around their common 
center of mass and with a negligible density of matter at 
the artificial container walls, as the one we describe e.g. 
in . On the other hand, if in our setting we release the 
constraint on the functions bi^m we get configurations 
similar to those they find. 

A possible way to enforce preservation of center of mass 
within molecular dynamics may be to take only those ini- 
tial conditions with positive parity, where for each parti- 
cle at {ri,Pi) there is a partner starting at opposite initial 
point in phase space (— r^,— p;). Such pairs of particles 
would hit the container simultaneously at opposite sides 
and the center of mass of the whole system without the 
container would not move. This should lead to symmet- 
ric binary configurations as the ones found by us in Qj. 
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